Exactly soluble model of resonant energy transfer between molecules 
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Forster's theory of resonant energy transfer (FRET) predicts the strength and range of exciton 
transport between separated molecules. We introduce an exactly soluble model for FRET which 
reproduces Forster's results as well as incorporating quantum coherence effects. As an application 
the model is used to analyze a system composed of quantum dots and the protein bacteriorhodopsin. 
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I. INTRODUCTION 

It has long been known that near-field electrodynamics 
allows energy transfer without emission of real photons. 
A striking example of this principle is provided by re- 
cent advances in wireless non radiative energy transfer 
1|; another recent example is Auger- mediated 'sticking' 
2| , whereby scattering states of positrons may transition 
to bound states in a metal, by transferring energy to a 
valence electron which can then leave the surface of the 
metal. Nowadays the phenomenon of non-radiative de- 
cay is of great scientific interest in many different fields of 
physics and chemistry [3]. In particular, new paradigms 
for solar energy conversion make use of nonradiative cou- 
pling for direct transfer of energy from the excitons cre- 
ated in the solar absorber to high mobility charge carri- 
ers Q , 01 • The mechanism for this transfer relies on the 
near-field resonance of electric dipoles, and is generally 
known as Forster resonance energy transfer (FRET) p- 
[1] . In its simplest formulation, FRET [9] is the quantum 
version of a classical resonance phenomenon, whereby os- 
cillating electric dipoles exchange energy through their 
mutual electric fields ^lOi] . Some studies in the quantum 
version have considered the possibility of coherent inter- 
actions between the dipoles [ll| - [l8| . Such quantum co- 
herence has been observed in the FMO complex, [l^ and 
it has been suggested that this may partly explain the 
high efficiency of energy transfer between chromophores. 
Interestingly, it has recently been demonstrated that the 
classical dipole model reproduces the quantum coherence 
as well [lo, [111- In this work, we revisit the analysis 
of the rate and efficiency of FRET, in the context of 
a donor and an acceptor species with comparable elec- 
tronic energy gaps. In this situation FRET is evidenced 
by decreased natural fluorescence from the donor, and 
enhanced fluorescence from the acceptor. The distance 
over which FRET has been observed ranges from 1 to 10 
nm, with the strength varying as the inverse sixth power 
of the separation. 

The fundamental mechanism underlying FRET is res- 
onance between excited electronic states in the donor and 



the acceptor molecules. The excited electronic states 
have non-zero electric dipole moments, and the result- 
ing dipoles experience a Coulomb interaction. The en- 
ergy exchange is complicated by the coupling of electronic 
states to vibronic molecular states, leading to a broaden- 
ing of the linewidths and a weakening of the resonant in- 
teraction. This vibronic coupling explains the difference 
between early inaccurate calculations of FRET efficiency 
(by Perrin [lOj and others) which were based solely on 
dipole resonance, and the later more successful calcula- 
tions by Forster Q which included vibronic effects. 

In this paper we consider a simple model for FRET 
which incorporates both electronic and vibronic effects. 
The model applies in situations where the donor molecule 
is rigid, with weak coupling between its electronic and its 
vibronic states, while the acceptor has strong electronic- 
vibronic coupling in its excited state. In this situation 
the model is exactly solvable and thus allows a compar- 
ison with the perturbative formulas derived by Forster 
and others. In particular we derive exact formulas for 
FRET efficiency and the Forster radius, and we com- 
pare these to the well-known Forster formulas. Further- 
more, the model is fully quantum mechanical and pre- 
dicts coherent oscillations between donor and acceptor 
under strong FRET conditions. The model contains a pa- 
rameter which determines the strength of the electronic- 
vibronic coupling in the acceptor, and for weak coupling 
the model reproduces the long-range interactions (up to 
100 nm) calculated by Perrin. [l^ 

As described above, the model applies to a FRET sys- 
tem where a rigid donor species, with weak coupling be- 
tween electronic and vibronic states, interacts with an 
acceptor species where electronic and vibronic states are 
strongly coupled. Thus the donor is modeled by a simple 
two-state system, corresponding to its electronic ground 
and excited states. The donor's vibronic degrees of free- 
dom are frozen and do not appear in the model. As we 
discuss in more detail below, this kind of model can be 
realized in practice with quantum dots (QDs). For the 
acceptor we again include only two electronic states, cor- 
responding to the ground and excited states, but in addi- 
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tion, wc include vibronic effects in the excited state. The 
excited band is described in detail below. For the mo- 
ment, we note that the electronic- vibronic coupling is de- 
rived from the Born-Oppenheimer approximation and as- 
sumes that the vibronic degrees of freedom are entrained 
to the electronic state. Initially we assume that the vi- 
bronic degrees of freedom are also frozen in the acceptor 
ground state. Later we indicate how nonzero tempera- 
ture effects may be included by unfreezing these degrees 
of freedom. 

The key step in the sohition of our model is the 
reduction to a finite dimensional system which ex- 
hibits the same efficiency for resonant energy transfer. 
The efficiency can be computed exactly for this finite- 
dimensional model, and this provides our exact results 
for the full model. In order to compare with Forster's 
formulas, we also compiite the absorption coefficient for 
our model, and we use this to evaluate the overlap inte- 
gral which appears in the standard rate formulas. 

As an application we apply our model to the analysis of 
one particular system of this type which exhibits FRET, 
namely the pairing of a QD donor species with bacteri- 
orhodopsin (bR) as the acceptor. The QD is known to 
have sharp emission lines for fluorescence and, thus, is 
a good candidate for the rigid donor molecule described 
above. We find good agreement with other calculations 
of the FRET rate for this system. 



II. DESCRIPTION OF THE MODEL 

We employ the Born-Oppenheimer approximation and 
assume that the electronic state of the molecule deter- 
mines its overall character, so that the vibronic state is 
entrained to the electronic state. Thus the state space 
is a direct sum Ho (B Hi (B ■ ■ ■ with one factor for each 
electronic state. The electronic configuration determines 
an effective Hamiltonian for the vibronic degrees of free- 
dom, and each space Hk is spanned by these vibronic 
states attached to the corresponding ground state. The 
spaces Hk may be discrete or continuous, depending on 
the striicture of the vibronic states. For simplicity we 
include only two electronic states, the ground state and 
the excited state, thus the state space is Hgr © Hexc- 



A. The ground subspace Hgr 

We assume a non-degenerate electronic ground state 
\ipgr)', and we also assume that the vibronic modes are 
frozen, so Hgr is one-dimensional. For the acceptor this 
assumes zero temperature. Later we extend to nonzero 
temperatures by including vibronic ground states. 



B. The excited subspace Hexc 

Again, we assume a non-degenerate electronic excited 
state \ijjexc)- For the donor, the vibronic modes are 
frozen, so that Hexc is one-dimensional. However, for 
the acceptor, the electronic state determines an effective 
Hamiltonian for the vibronic states, which are labeled 
by their energy eigenvalues e. For this subspace we as- 
sume that the vibronic states form a continuous band 
with a uniform density of states, with eigenvalues extend- 
ing from — (X) to -|-oo. Under this assumption we ignore 
any edge effects in the band. Thus the space Hexc is iso- 
morphic to the one-particle Hilbert space -L^(R), and we 
represent a state as a square integrable function (jjexds) 
where 

/oo 
\(l>exc{e)\^de=l . (1) 
-CXD 

The time evolution is (j)exc{£) — ^ e~*'^*0exc(e)- For conve- 
nience we denote the Hamiltonian in this basis as h, so 
that 

{h(f)exc){s) =S(t>exc{s) ■ (2) 

C. Excitons 

As long as the electronic state does not change, the 

dynamics of the vibronic state is completely determined 
by the fixed effective Hamiltonian corresponding to this 
electronic configuration (here we are neglecting any feed- 
back reaction from the vibronic modes on the electronic 
modes). However, when the molecule undergoes an elec- 
tronic transition, for example by photon absorption, the 
effective Hamiltonian for the vibronic states immediately 
changes. This sudden change creates an excited vibronic 
state, as the previously stationary vibronic state becomes 
a superposition of energy eigenstates states of the new 
Hamiltonian. We call this vibronic state an exciton. The 
exciton behaves like a delocalized one-particle state. In 
our model, exciton states will arise only in the excited 
band of the acceptor, due to a transition from the ground 
state. We assume an average energy £^3 for these exciton 
states. 



D. Transitions 

Turning now to transitions, we consider only radiative 
interactions which act solely on the electronic state. Thus 
transitions of the vibronic state occur as a consequence 
of the change of the effective Hamiltonian due to the 
electronic transition. The electronic matrix element due 
to the interaction V is 

{i^exc\V\ijgr) . (3) 

We find an explicit form for this matrix element for the 
situations of interest, namely, direct photon absorption 
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and resonant excitation through the Coulomb interac- 
tion. In order to determine the vibronic matrix element, 
we follow Jortner and proceed by analogy with the 
derivation of the lineshape of a resonance. Recall that a 
resonance is a perturbation of an embedded eigenvalue in 
continuous spectrum. The perturbation causes the eigen- 
value to dissolve, accompanied by the emission of a one- 
particle state in the continuous spectrum. In our model 
this one-particle state is the exciton. Thus the transi- 
tion from a vibronic ground state is accompanied by the 
creation of an exciton, which is a normalized excited vi- 
bronic state. The Breit-Wigner form for the lineshape 
of the resonance is a Lorentzian, where the width cor- 
responds to the lifetime of the resonance, [l^l and the 
center is the average energy. We assume the same form 
for the exciton, so the wave function of the exciton (in 
the diagonal energy representation) is the square root of 
a Lorentzian: 



electronic parts of the states \D*A) and \DA*). Thus the 
electronic transition matrix element is 



fie) 



27r (£ - ^3) 



+ f7 



(4) 



Using a Lorentzian form implicitly assumes that the en- 
ergy band extends from — cxd to +00. We make this as- 
sumption, thus ignoring any edge effects in the band. The 
width 7 depends on the particular system and determines 
the lifetime of the exciton. The Lorentzian is centered at 
energy E3, corresponding to the average exciton energy. 
We also include a phase factor e'^ , which may depend on 
e. As we will see this phase factor is irrelevant to the 
calculation of the FRET efficiency. 



III. DYNAMICS OF THE MODEL 

The FRET interaction between the donor and the ac- 
ceptor causes an exchange of energy as the excited state 
is transferred from one to the other. Although the trans- 
fer becomes irreversible after some time, the initial in- 
teraction is unitary and thus admits the possibility of 
oscillations between donor and acceptor. 

We use standard notation, with D for donor ground 
state, D* for donor excited state, A for acceptor ground 
state, and A* for acceptor excited state. Thus the cou- 
pled donor-acceptor system is described by the collection 
of states 



IDA), \D*A), \DA*), \D*A*), 



(5) 



where both \DA) and \D*A) are single states, while 
\DA*) and \D*A*) contain the exciton subspace. 



A. The Hamiltonian 

The FRET interaction is dipole-dipole in lowest order, 
and so its strength decays as the inverse third power of 
the distance. The Hamiltonian is determined by the ma- 
trix element of the Coulomb interaction Vc between the 



U = {DA*\Vc\D*A) 

= ±fr>D-BA- ^{-Dd ■ R)(D^ • R; 



(6) 



where R is the separation between the systems, and 
and Da are the transition dipole moments of the donor 
and acceptor, respectively. We use atomic units through- 
out, and we assume that the dielectric constant is 1. We 
look in detail at specific models later, but for the mo- 
ment we note that for typical systems the dipole moment 
is about 10 D, so at separations of around 5 nm the in- 
teraction energy U 10~^ eV, compared to the typical 
energy gaps between ground and excited states of 1 — 2 
eV. 

In the absence of other effects, we could analyze the 
dynamics of this coupled system by restricting to the 
subspace spanned by the states \D*A) and \DA*), and 
computing the time evolution of the initial state \D*A) 
under the influence of the interaction U. The Hamilto- 
nian is 



H 



El 
U\f) 



U{f\ 
E2 + h 



(7) 



where Ei and E2 are the electronic energy gaps of the 
donor and acceptor, respectively; U is the interaction 
matrix element defined in Eq. (|6|); and h is the diagonal 
energy operator of the continuous exciton band in the 
excited state. Also, |/) denotes the creation operator for 
the exciton as in Eq. (|4]), and (/| denotes the correspond- 
ing annihilation operator. 



B. The master equation 

In our model we also include the effects of fluores- 
cent decay from the excited state to the ground state 
for both systems. We do this by introducing jump op- 
erators for the (irreversible) fluorescent decays from ex- 
cited to ground state and use a master equation to com- 
pute the time evolution of the density matrix. So we 
are using the Markov approximation for the couplin g to 
the electromagnetic field which causes fluorescence. |23| 
In order to separate the outcomes from the two excited 
states \D*A) and \DA*), we use two copies of the ground 
state to indicate which molecule has decayed (since the 
jump operators are irreversible, there is no coupling from 
these ground states back to the excited states). Thus 
the system is represented by a density matrix p, span- 
ning the states \D*A), \DA*), \DA), and \DA), where 
\DA) and \DA) are copies of the ground state, and the 
two fluorescent decay channels are \D*A) \DA) and 
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\DA*) \DA). In this subspace the Hamiltonian is 



IV. SOLUTION OF THE MASTER EQUATION 



H 



DA 



f E, U(f\ 0\ 

U\f) E2 + h 



\ 0/ 



(8) 



The effects of fluorescence are implemented by the jump 
operators: 



Ji = 



/O 0\ 














01 














1 0' 












Vo 0/ 






{Ui\ 





0/ 



(9) 



where form an orthonormal basis in the exciton 

space. The master equation is 



^ = -i[HDA,p] 



(10) 



71 
' 2 



(2Jipjl-jlJip^pjlJi 



where 71, 72 are the rates for fluorescence \D*A) \DA) 
and \DA*) \DA), respectively. 



C. Definition of tlie efficiency, Forster radius and 
FRET rate 



Using the master equation in Eq. (|T0|) . it is possible 
to compute the probability of exciton transfer from the 
donor to the receiver and to find the efficiency of this 
process. In the absence of FRET, the system ultimately 
ends up in state \DA). Thus we define the efficiency of 
FRET to be the long-run probability that this does not 
happen; that is, 



F ^1- lim {DA\p{t)\DA) 

t^OC 



(11) 



The Forster radius Rq is then defined by the condition 
that at this separation the efficiency reaches 50%. That 
is, 



Ro = max{i? : F > 0.5}. 



(12) 



The FRET rate 'Jfret can also be computed from the 
efficiency, by comparing the rates for fiuorescence of the 
donor and FRET: 



IFRET = 71 



1 - F 



(13) 



where 71 is the natural fluorescence rate for the donor. 
This follows from the relation 



F 



IFRET 

If RET + 71 



(14) 



We solve the master equation given in Eq. (|T0)) with 
the initial condition p(0) = \D* A){D* A\, corresponding 
to the donor in its excited state and the acceptor in the 
ground state. The solution of Eq. (|T0| takes the block 
diagonal form 



Pit) 

The top left 2x2 block is 



Pii(t) Pi2{t) 

P2l(<) P22{t) 

P33(i) 

Piiit)/ 



(15) 



Pll{t) Pl2{t) 
P2l{t) P22{t) 



(16) 



where B is the non-Hermitian operator acting on the 
space spanned by \D*A) and \DA*{e)), i.e.. 



B = 



El- ^11 U{f\ 

U\f) E2+h-^^2j' 



(17) 



and I/) is the Lorentzian function in Eq. Q. The re- 
maining diagonal entries in Eq. (1151) are given by 



Psiit) = 71 / Piiis)ds 
Jo 

Piiit) = 72 / Trp22(s)ds. 
Jo 



(18) 
(19) 



In order to facilitate the notation, define the state 

iV-o) = (J) . (20) 

It follows that the efficiency is given by 

F=l-p33(oo) = l-7i / |(^o|e-^^1^o)|' rfs.(21) 

"'0 

Thus the calculation reduces to the problem of finding 
matrix elements of the operator e^'^*. This is straight- 
forward because i? is a rank 1 perturbation of a diagonal 
operator. The key step is the reduction to a related two- 
state system. Namely, define the 2x2 matrix 



B = 



i?i-f7i 
U E2 



U 



^3- 5(7 + 72) 

It is shown in Appendix A that 
and thus 



F = l--fi 



ds. 



(22) 



(23) 



(24) 



Thus the effect of the exciton coupling in this model is 
the same as in a two-state model with a second channel 
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for decay of the excited state to the ground state, at a 
rate which is the inverse hfetime of the exciton. 

The efficiency F in ((24| can be evaluated by finding 
the eigenvectors and eigenvalues of B. The assumption 
that 71 , 72 > implies that B has two eigenvalues with 
negative imaginary parts. These can be computed explic- 
itly, and the initial state j-i/'o) can be written as a linear 
combination of the eigenvectors. The result is 



(l + r)|C/|^ 



where 



and 



(1 + r)2|C/|2 + 4r [Ei + - E2 - E^y] 



r = 71 + 72 + 7 



71 



72 + 7 



(25) 



(26) 



(27) 



Forster radius and FRET rate 



From Eq. (|25p we compute the Forster radius and the 
FRET efhciency for this model. Setting F = 1/2 we get 
the condition 



4r 



71+72+7 



(28) 

{El — E2 — E^Y 

. Introducing an 



The interaction U is given by Eq. ^ 
angular factor n this gives 

and, hence, the formula for the Forster radius 



4r 



16 1 7i + 72 + 7 



(29) 



■(30) 



{El -E2- E^f 



B. Exciton lifetime and scaling at resonance 

At resonance where the energies match we have Ei 
E2 — E3 — 0. In this case we have 



F 



72+7 



\U\' 



71+72+7 |C/|2 + 2lh^ 
Define the dimensionless parameter 

_ 4|[/|2 
^ 71(72+7)' 



(31) 



(32) 



then at resonance we get 



F = 



72+7 



71 + 72 + 7 '/ + 1 



(33) 



It is reasonable to expect that the inverse exciton lifetime 
7 is much larger than the fluorescence rates 71 , 72 . In this 
case the FRET efficiency at resonance takes the simple 
form 



F = 



V 



4|t/|2 
71 7 



(34) 



This same scaling relation has been found for resonant en- 
ergy transfer between classical oscillators [H, and seems 
to be a general feature of this type of phenomenon. 



C. Temperature dependence 

In order to incorporate non-zero temperature effects, 
we introduce vibronic ground states for the acceptor, la- 
beled {\4'gr{£i))}, where ei is the energy in the ground 
state vibronic band. These states are assumed non degen- 
erate, and thus completely labeled by their energy eigen- 
value El, with respect to the Hamiltonian determined by 
the electronic state {ipgr)- Thus a general state in the 
acceptor band will be 



|A)^^c(ei)|0,.(ei)), 



(35) 



with the normalization condition |c(£i)|2 = 1. The 
vibronic Hamiltonian is diagonal in this representation, 
and so the time evolution of a state is 



lA)^^c(ei) 



4>gr{ei)). 



(36) 



To include nonzero temperatures, we assume a Boltz- 
mann distribution for the initial equilibrium state; that 
is, 



Pe 



<t>gr{£l)){(t>gr{£l)\ 



(37) 



We then compute the thermal average of the efficiency 
over initial vibronic energies. Thus the temperature- 
dependent efficiency is 



F{T) 



F{ei)dei 



(38) 



where F{ei) is given by Eq. (|25|) with Ei replaced by the 
initial energy Ei+Si. By using this formula for F{ei) we 
assume that initially the acceptor is in a pure vibronic 
ground state \4>gr{£i)), which transitions to the exciton 
state /(e) due to the FRET interaction. Our main sim- 
plification is the following: we assume that the reverse 
operation causes a transition from the exciton state back 
to the same initial vibronic state \(j)gr{£i)). So under 
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this reverse operation an exciton state ijj{e) is mapped to 
a\(pgr{si)), where the amphtude is a = J f{e)'ip(e2)de. 
Thus the acceptor always returns to its initial ground 
state (this assumption has also been for exact calcula- 
tions of coherent exciton scattering fiij ) . 



Carrying out the summation in Eq. (138)) requires 
knowledge of the phonon spectrum of the acceptor, and 
we do not pursue the question further here. However, we 
note that this effect of the temperature is expected to be 
small because at room temperature kgT is significantly 
smaller than the energy scale given by 7 ^ 0.1 eV. 



where is the Bohr radius. Using this expression for 
the absorption coefficient in Eq. (15^ and setting lu = Ei, 
we get (in a.u.) 



9k^Q 



D 



C 

El 



(44) 



{Ei-E2-E3)^+jyA 



El 

El 



Turning now to our Eq. (j30p . we use the Einstein A 
coefficient to relate the fluorescence rate and the dipole 
strength of the donor: 



V. COMPARISON WITH STANDARD FRET 
FORMULAS 

The standard Forster formulas for FRET rate and 
Forster radius involve the overlap integral between the 
normalized donor fluorescence spectrum and the accep- 
tor absorption spectrum. The formula for the Forster 
radius (in cm) is [1| 



6 _ 90001n(10)K2QD , 



(39) 



where is the donor quantum yield, Na is Avogadro's 
number, and J is the overlap integral (in cm^dm'^/mol). 



J 



(40) 



(recall that we have assumed that the refractive index 
is 1). Here Fd(A) is the donor emission spectrum [nor- 
malized so that Fd{X) dX = 1, where A is the wave 
length (in cm)], and ca is the molar absorption coefficient 
in (cm~^ dm^/mol). Inserting values for the constants 
gives 



R% = 8.79 X 10-25 Qd J- 



(41) 



In our two-level model the donor is assumed to have a 
narrow band fluorescence spectrum centered at energy 
El, so the overlap integral is essentially caIX) X'^ evalu- 
ated at the wavenumber A = 2ttc/Ei. Thus the formula 
is 

= 8.79 X 10^25 ^2 ^^(^2nc/Ei) {2nc/Ei)-^. (42) 

The acceptor's molar absorption coefficient can be 
evaluated from the standard absorption rate for the tran- 
sition from ground state to excited state, using Fermi's 
Golden Rule to compute the rate. The details are carried 
out in Appendix B and the result (in cm"^ dni"^/mol= 10 
m^/mol) is 



eAln(lO) 27r 



10ElD\^ 



ujc {uj~E2-E3)^+-f^/4' 



(43) 



71 



3 c3 



(45) 



This approximation for 71 given by Eq. (|45p neglects 
spectral shifts and broadening. Thus our expression for 
Rq becomes 



(72 + 7)D\ 



(46) 



{El -E2- E^f + (71 + 72 + 7)V16 ■ 

Comparing Eq. (|46p and Eq. pij) . we see that our re- 
sult differs from the standard one in several ways. In 
particular, the width of the Lorentzian is (71 -I- 72 -f 7)/4 
instead of 7/2, and the factor {E2/Eif' is absent in 
Eq. However, for realistic models we expect that 

7 >> 71,72, hence r << 1. So when we compare the 
values at the resonant energy (where Ei — E2 — E^ = 0) 
and if we set = 1, then we get the ratio 



Rq 

R^ 



= 0.7783 



El 
E2 



1/3 



(47) 



We expect that E^ ~ £'2/10 (the energy difference be- 
tween the donor's emission peak and the acceptor's ab- 
sorption peak), and thus E1/E2 ~ 1.11, so we find a quite 
close agreement with the standard result. 



VI. APPLICATION: QUANTUM DOTS AND 
BACTERIORHODOPSIN 

Recent proposals for improved dye-sensitized solar cells 
(^S*, '26] involve replacing the liquid dye by nanoparticles 
attached to a substrate and exploiting FRET to achieve 
efficient energy transfer . One candidate material is a 
mixture of QDs and the protein bacteriorhodopsin (bR) 
p8l |. [2^, [201 • In one scenario the QD would act as an 
antenna for photon absorption, with subsequent transfer 
to the retinal complex in bR. The retinal complex in bR 
is known to be an efficient absorber of photons through 
direct capture, and this same efficiency is expected for 
non-radiative transfer of excitons from QD to bR via the 
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FIG. 1: (Color online) Schematics of the continuum model 
for FRET. The donor is a quantum dot while the acceptor is 
an opto-electric protein such as bacteriorhodopsin. 



FRET mechanism. Q The methods developed in this pa- 
per can be used to evaluate the efhciency of FRET in 
this hybrid system. The QD has a band gapof approx- 
imately 2 eV (depending on its diameter Vj3lL l32l | and 
after photon absorption it rapidly relaxes to its lowest 
energy excited state, thus the QD is well modeled as a 
two-state system. 

In its ground state the retinal molecule has a planar 
conformation. Upon excitation it briefly enters a band 
of planar excited states (due to an electronic transition 
consistent with the Franck-Condon principle), and then 
rapidly relaxes to a non planar conformation. [33j The lat- 
ter transition occurs within a few hundred femtoseconds, 
is effectively irreversible, and thus signals the transfer of 
the excitation to bR. The planar excited state lies in a 
band of closely spaced levels corresponding to different 
vibrational and rotational states. Recent studies have 
demonstrated that coherence persists in the exciton state 
for several hundred femtoseconds after initial excitation 

Our model QD/bR is schematized in Fig. 1 and it is 
a simplified version of the more general model described 
above. The values of the various parameters can be ob- 
tained from known properties of the molecules. The wave 
function / has the Lorentzian form, centered at the ex- 
citon energy E3. We assume transfer on resonance, so 
that El = E2 + E3 = 2 eV. The width 7 is the inverse 
lifetime of the exciton, which is known from coherence 
analysis to be at least 100 fs, so 7 is upper bounded by 
around 0.05 eV. The rate 71 is set by the QD experimen- 
tal lifetime, [34] so 7f ^ — 16 ns, and the rate — 500 
fs.[33j However, they are not important in the calcula- 
tions since they are much smaller than 7. The FRET 
coupling strength U is determined by the formula 16] 



U = 



K Dqd DbR 



(48) 




0.5 





Y-0.05 eV 
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FIG. 2: (Color online) (a) Efficiency for QD-bR as a func- 
tion of separation; (b) occupation probability of initial excited 
state as a function of time. Arrows show the position of Ro 
for two values of 7. 



dipole moments of the QD and bR respectively, and the 
angular factor k depends on the orientations of dipoles 
relative to the separation between molecules. We use 
values Dqo — Dhu = 10 D, — 1 (permittivity of 
medium, assumed dry), and k = 1 and keep the separa- 
tion distance i? as a free parameter. In atomic units this 
gives 



U = 



15.479 



i?3 



(49) 



where is the permittivity of the medium, R is the sep- 
aration between the molecules, Dqo and Dbn are the 



Figure 2 (a) shows the efficiency as a function of R for 
these values. The FRET distance Rq is consistent with 
the one estimated from experiments of about 7 — 8 nm 
[soj . The curve almost exactly matches the phenomeno- 
logical formula [^-[8] for efficiency E = r{/{R^ + R^). 
Figure 2 (b) shows the occupation probability of the ini- 
tial donor excited state as a function of time, for the 
same parameter values and with a separation R = 2 nm. 
Coherent oscillations are apparent if 7 = 0.05 eV. 



VII. CONCLUSIONS 

We have introduced an exactly solvable model for 
FRET, which captures the key features of Forster's 
electronic-excitonic coupling in a microscopic quantum 
mechanical setting. The standard Forster equations are 
accurate when the following conditions are satisfied: 

a) The dipole-dipole approximation for the electronic 
coupling can be employed appropriately for the 
donor-acceptor interaction. 

b) Interactions among donors or acceptors and static 
disorder effects leading to spectral line broadening 
can be neglected. 

c) The energy transfer dynamics is incoherent. 
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Our approach goes at least beyond condition c since, at 
short distances and times, the formaUsm is able to cap- 
ture the coherent energy transfer. The model uses a mas- 
ter equation with Lindblad operators to take account of 
fluorescence and relaxation effects, and uses a continuum 
of excited states in the acceptor to implement the exci- 
ton dynamics. The model is robust and can easily be 
extended to include more complicated exciton dynamics. 
As a concrete application the model is used to analyze 
FRET coupling between a QD and bR, where it makes 
realistic predictions of the FRET distance. 
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Appendix A 

We use resolvent techniques to compute the exponen- 
tial of B. Recall the resolvent representation 

1 
27ri 



-iBs 



{z-By'^dz, 



(Al) 



where the line integral encloses the spectrum of B in the 
complex plane. The assumption that 71,72 > implies 
that B has no spectrum in the upper half-plane, so the 
resolvent (z — B)~^ is analytic in the upper half-plane. 
There is a cut along the real axis where the spectrum of 
/i2 lies, and also, possibly, poles in the lower half-plane. 
The line integral encloses the cut and also any poles in the 
lower half-plane. In the lower half-plane the integration 
contour can be deformed to a large semicircle z — Re~^^ 
with < 9 < TT. For s > the contribution of this 
semicircle vanishes in the limit i? — > cxd, thus for s > 
the line integral in Eq. (|A1[) can be written as 



1 

1 

27ri 



oo+te 



e-^''" [z ~ By^ dz (A2) 



— oo-\~ie 



1 /"oo-t-ie 

— / e-"" {B - z)-^ dz. 

J-oo+ie 



This leads to the formula 

-1 />oo+?e 

(V'ole-^^^IVo) = ^ / (A3) 



(^o|(S-z)-Vo)d^. 



We next derive an explicit formula for the matrix ele- 
ment {tpoKB — z)~^\'ipo) appearing on the right-hand side 
above, under the assumption that Imz > 0. Define 



yhiiz) h2{z)) 
Then the Feshbach method yields 

(^o|(S-2)"Vo) = 



where 



= \E,---n-z-\UYM 



M^{f\{E2 + h--^2-z)-^\f) 



(A4) 



(A5) 

-1 



(A6) 



Using the Lorentzian form for /, and the diagonal energy 
operator ft,, the matrix element M is 



M = — 
27r 



1 



, fc - ^3 + §7 
1 1 



(A7) 



fc - S3 - |7 ^2 



■dk. 



■72 - z 



This integral may be computed by completing the con- 
tour in the lower half-plane and evaluating the sum of 
the residues. For Imz > this gives 



M 



1 



E2 + Ez- ^y + -12) - z 

Inserting this into Eq. (|A5P leads to the expression 



(AS) 



/ii(z) = (£'1 - -71 - z 



(A9) 



£'2+£'3-5(7 + 72)-z 

The key observation now is that Eq. (IA9p is the resolvent 

of the reduced two-state system defined by the matrix 
introduced in Eq. (f^ : that is, 

'^1 - f 71 U 



B = 



(AlO) 



U -B2 + S3 - 5 (7 + 72), 
It follows that 

{^ol(i? - z)-^\^^) = {M{B - z)-Vo) , (All) 
and hence we obtain for all s 

(V'ole-^^IVo) = (V'ole-'^IV'o) ■ (A12) 

Appendix B 

The acceptor's molar absorption coefficient can be 
computed by using the following time-dependent Hamil- 
tonian for the electronic transition in the presence of a 
classical field: 



H{t) - ( Jl 



ve 



(Bl) 
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Here E2 is the energy of the excited state, uj is the fre- 
quency of the radiation, and z/ is the couphng between 
the field and the system. This coupling is given by 



ie f 
mc J 



(B2) 



where 4'exc and V'gr are the excited and ground state 
wave functions, and Aq is the field strength (assumed 
constant). We have used the rotating wave approxima- 
tion and dropped the coimtcr-rotating term proportional 
to e*"^*. Standard dipole approximations lead to 



!/= \A0\E2DA . 

c 



(B3) 



The radiation intensity (power per unit area) is related 
to the field strength through the time-averaged Poynting 
vector, and this gives 



OJ 



2nc 



(B4) 



Combining the electronic transition rate with the exciton 
amplitude, Fermi's Golden Rule gives the following rate 
for transitions in this radiation field: 



The absorption coefficient determines the rate of energy 
absorption by an ensemble of molecules. Consider a slab 
of absorber with unit cross-sectional area, and thickness 
X, illuminated by light of frequency oj and intensity 
Then the output intensity is given by the Beer-Lambert 
law 



(B6) 



where ca is the molar absorption coefficient and k is the 
concentration (in moles per unit volume) . The number of 
molecules in the slab is N^kx, where A''^ is Avogadro's 
number. The energy absorbed per unit time by tran- 
sitions is thus NAkxaco. Equating this to the energy 
difference between input and output gives the result in 
atomic units: 



e^ln(lO) = 



Na<7C0 



NACo2Tr\iy\^\f{cu-E2)\' 



NaEId\j 



2tt 



wc {UJ-E2- Esf + i^l^ 



a = 2'K\vf\f{u-E2t 



(B5) 
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